T c suppression in co-doped striped cuprates 
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We propose a model that explains the reduction of T c due to the pinning of stripes by planar impurity 
co-doping in cuprates. A geometrical argument about the planar fraction of carriers affected by stripe 
pinning leads to a a linear T c suppression as a function of impurity concentration z. The critical 
value 2 C for the vanishing of superconductivity is shown to scale like T% in the under-doped regime 
and becomes universal in the optimally- and over-doped regimes. Our theory agrees very well with 
the experimental data in single- and bi-layer cuprates co-doped with Zn, Li, Co, etc... 

PACS numbers: 74.20.Mn, 74.50. +r, 74.72.Dn, 74.80.Bj 



One of the most striking properties of high temper- 
ature superconductors (HTSC) is the sensitivity of the 
critical temperature, T c , to planar impurities which are 
introduced with the substitution of the Cu atoms in the 
Cu02 planes. T c is suppressed with a few percent of dop- 
ing almost independently of the magnetic nature of the 
impurity [Q . In particular it has been shown experimen- 
tally that the HTSC undergo a superconductor to insu- 
lator transition due to Zn doping |g{ . The fast depression 
of T c with Zn has been studied extensively by NMR, /xSR 
HQ , and infrared techniques || and a heated debate has 
been developed on the interpretation of the experimental 
data It has also been established that this suppres- 
sion is more robust in the under-doped compounds 0. 
The reduction of T c has been assigned to formation of 
magnetic defects [|]||, electron scattering by disorder in 
the presence of a d-wave order parameter p~0| . |TT| ] , unitary 
scattering jl2| , local variation of the superconducting gap 
fL3| |, and weak localization jwj. However, although the 
existing theories describe well the behavior of these sys- 
tems at low impurity concentration, they systematically 
deviate from the experimental data at higher doping con- 
centration. 

In this work we introduce a scenario for the destruc- 
tion of superconductivity which is based on the stripe 
picture. We propose that Zn pins stripes and slows 
down their dynamics in a region of size R in its sur- 
rounding. This effect resembles the "swiss cheese" -like 
model with Zn sites creating voids in the superconduct- 
ing mesh and destroying superfluid density In our 
model, however, the superfluid density remains unaltered 
upon Zn co-doping, whereas the stripe inertia increases 
locally. Since the stripe mass density p is inversely pro- 
portional to the hopping energy t |l5|], a local decrease 
of t by pinning will enhance the stripe effective mass, p. 
Magnetic studies of the effect of Zn doping in cuprates 
indicate that co-doping indeed leads to charge carrier 
localization p 16 171]. T c {x,z) is a function of the car- 
rier density x and z and it is determined by the ratio 



between the 3D superfluid density n s and the a6-plane 
charge carriers effective mass m* b . Within the stripe 
picture we must have T c oc n s /p. Hence, the pinning 
of stripes in the vicinity of the impurity implies a re- 
duction of T c by the enhancement of p. As a result 
of our model we obtain T c (x, z)/T c (x, 0) = 1 — z/z c (x), 
where z c (x) is the critical concentration required to com- 
pletely suppress superconductivity. In under-doped com- 
pounds, z c (x) oc x 2 and scales with T c 2 (cc,0) in con- 
trast with the case of ordinary superconductors. In 
optimally- and over-doped materials z c (x) is indepen- 
dent of x and exhibits universal behavior. Our results 
give a very good description of the experimental data 
in La 2 - :E Sr :E Cu04 (LSCO), YBa 2 Cu 3 7 _ e (YBCO), and 
Bi 2 Sr 2 Cai_ y Y y Cu208+e (BSCO) co-doped with different 
impurities, such as Zn, Li, Co, etc. 
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FIG. 1. Effect of Zn on the stripe grid. Dark dashed lines 
mark the regions of stripes with hnite superfluid density. The 
circle shows the region where hopping is suppressed and the 
superfluid effective mass is increased. 

Striped phases of holes have been observed experimen- 
tally in the superconductor Lai.6-xNdo.4Sr x Cu04 |Tg| ] 
and in La 2 _ 2 ,Sr :E Ni04+ a | fi9| , which is an antiferromag- 
netic insulator. Moreover, stripe formation provides a 
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simple explanation for the observed magnetic incommen- 
surability in LSCO §0j and YBCO §|. The magnetic 
incommensurability appears in neutron scattering by the 
splitting of the commensurate peak at Q = (Tt/a,ir/a) 
(a m 3.8A is the lattice spacing) by a quantity S, which 
is inversely proportional to the average inter-stripe dis- 
tance £ p^ ]. In under-doped LSCO, the charge stripes 
behave as an incompressible quantum fluid and S is pro- 
portional to x p3|. In this regime T c scales with 5, p2] 



T c (x) 
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where £ is given in A. On the other hand, in the opti- 
mally and over-doped compounds the incommensurabil- 
ity 6 saturates to a constant value and £ becomes nearly 
independent of x |22fl . To date, there is strong evidence 
for stripe formation in LSCO compounds |22|Jlq| . Re- 
cent ion-channeling p4], NMR /xSR |26| and neu- 
tron scattering data fl27f also indicate that stripes might 
be present in YBCO, whereas a lack of experimental data 
still persists for BSCO. In what follows we will concen- 
trate on the physics of HTSC co-doped with planar im- 
purities, like La2-j;Sr x Cui_ 2 Zn z 04. We will verify that 
our results can describe equally well YBCO and BSCO 
co-doped with Zn, Li, and Co. 

For simplicity, we assume that the Zn atoms are lo- 
cated half distance between superconducting stripes and 
suppress stripe fluctuations within a circle of radius R 
around their position, as shown in Fig. [l]. This distance R 
is assumed to scale with the average inter-stripe distance 
£. Thus, we parameterize R = 7^/2 with the doping in- 
dependent phenomenological parameter 7 > 1. As it is 
shown in Fi g, p], the st ripe length which is pinned per Zn 
atom is A^/R 2 — £ 2 /A = 2£y/j 2 - 1. If we assume that 
all the Nzn Zn atoms take part in pinning the stripes, 
the total pinned length is Nz n 2£\/^ 2 — L The trans- 
verse kinetic energy density of the stripe is t/a, where t 
is the single hole kinetic energy. Thus, the suppressed en- 
ergy density in the plane is St = (t/a)2£y^j 2 — INzn/L 2 
where L is the sample size. Defining the impurity fraction 
z = Nz n {a/ L) 2 and the planar kinetic energy density 
t = t/(a£), we find that the suppressed energy density 
reads 



The energy density of the co-doped system is given by 
t(x, z) — t(x, 0) — St, which leads to 



t(x,z) _ z 



t(x,0) z c (x) 



where 



z c (x) 



(3) 



(4) 



Superconducting long range order, and therefore T c , is 
obtained when the stripe array attains phase coherence. 
This coherence can be achieved by the Josephson cou- 
pling between stripes which may occur via exchange of 
Cooper pairs |2^] or d x 2_ y 2 bosons |2j|. In this case the 
finite transition temperature is in the 2D XY universality 
class (stabilized by the inter-planar coupling). Thus, T c 
is given by 



2k B m* ab 



(5) 



where L c is the inter-planar distance. Notice that n s ~ 
l/(2a£(x)L c ) (Ref. (l8)) and m* ab is proportional to the 
stripe linear mass density p = h 2 /(ta 3 ) |]|. Therefore 
we conclude that 



rp t \ t(z) . 

T c {x,z) oc — — cx t(x,z). 
£{x) 

Combining (^|) with (^) we find: 

T c (x,z) t(x,z) z 



T c (x,0) t(x,0) z c (x) 



(6) 



(7) 
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FIG. 2. T c (x, z)/T c (x, 0) versus the Zn doping concentra- 
tion z for the lanthanate-compound. Dashed and dotted lines 
correspond to Eq. ((?]) with z c given by Eqs. (^) and (^|), re- 
spectively. Red circles denote x= 0.10 [BOj and the opened 
symbols are: + x = 0.15 (3(|, circles x = 0.15 |l|, lozenges x 
= 0.18 @, triangles x = 0.20 Q, squares x = 0.20 |o|. 

It is very illuminating to compare the expressions for 
z c (x), that is, the critical Zn doping for which T c van- 
ishes, in the various regimes. For under-doped lan- 
thanates (x < 1/8), neutron scattering data show that 
T c (x,0) oc 6 ||. Comparing (g) with ([[]), one finds an 
unforeseen result, which is 



z c (x) 



a(A)V T 2 (x,0) 
805 J T 



(8) 



in the under-doped regime. Moreover, in the optimally 
and over-doped compounds the average stripe separation 



2 



saturates to a value £ s , and therefore we expect z c to be 
an universal constant, independent of x, 
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FIG. 3. z c x T c : the continuous line is given by Eq. (^). 
The data points are experimental values for under-doped 
La-compound: black squares (Ref. (ll]]); blue circles (Ref. 
||), red circle (Ref. §). 

In Fig. [2], we plot the experimental data for LSCO co- 
doped with Zn J|^] . The dashed and dotted lines are the 
theoretical result (R) with z c given by Eqs. (||) and (||), 
respectively. The agreement between theory and exper- 
iment is excellent. One can easily observe that indeed 
the data for different x collapse into a single line for op- 
timally and over-doped compounds. In addition, in Fig. 
|3| we show the experimental data for z c x T c , as well as 
the theoretical prediction In both figures, we used 
7 ~ 1.42 so that R ~ 0.71£. Hence, for the optimally- 
and over-doped materials, for which the mean stripe dis- 
tance £ s = 4a = 15. 2A, we estimate R ~ 10A. The value 
of R can be independently obtained in neutron scattering 
experiments in Zn doped samples. At this point in time, 
however, we are not aware of any data in this regard for 
LSCO. 

The excellent agreement between theory and exper- 
iment is however not restricted to the Zn doped lan- 
thanates. In Fig. ^ we display the experimental data 
for YBCO (red points) co-doped with Li (squares) and 
Zn (circles) |33|j, a s well as for BSCO co-doped with Co 
(black points)p2| . The universal behavior at high dop- 
ing is very clear (Fig. |]a) as well as the different slopes 
characteristic of the under-doped compounds (Fig. ^b). 
Observe that YBCO and BSCO exhibit exactly the same 
features in terms of Zn doping, with the same value of 
z c in the over-doped regime. We assign this behavior 
to the fact that both systems are bi-layers while LSCO 
is a single-layer system. Moreover, it has been shown 
recently that the relationship of T c and 8 for YBCO is 
similar to the one in LSCO p3 . Using the value of the 



incommensurability 6 = 0.11 measured by neutron scat- 
tering in YBCO6.6 p7fl , which has T c = 60K, we can de- 
termine the stripe separation in this under-doped com- 
pound, £{x) ~ 4.55a. Then, by fitting the data from 
Ref. for the same compound (red circles in Fig. 4b, 
z c = 0.04) to Eqs. (0) and (|), we determine the only 
free parameter in our model, 7 — 1.17, corresponding to 
R ~ 10 A. Now, we use this value of 7 to analyze the over- 
doped regime. The universal line in Fig. 4a defines the 
critical doping in over-doped bi-layer cuprates, z c ~ 0.13. 
By replacing z c and 7 in Eq. , we can estimate the sat- 
uration value of the average stripe distance. We then find 
£ s ~ 3a, corresponding to R = r y£ s /2 ~ 6 A. These re- 
sults agree with magnetic measurements in over-doped 
YBCO with Zn impurities: NMR studies suggest the ex- 
istence of two different relaxation times ascribed to 63 Cu 
sites away from and near Zn impurities S and inelastic 
neutron scattering data J35| indicate that Zn induces a 
magnetic perturbation on a range of R ~ 7 A, very close 
to our estimates. Moreover, we predict that the satu- 
ration of the incommensurability in YBCO should start 
at optimal doping, where S s ~ 0.2 and £ s ~ 3a, in con- 
trast with the case of LSCO where the saturation starts 
at x ~ 0.12 and £ s ~ 4a. This result agrees with recent 
neutron scattering data in YBCO |34|] and establishes a 
very important relationship between two different kinds 
of experiments, namely, magnetization and neutron scat- 
tering. By determining the doping concentration x be- 
yond which the T c (x, z)/T c (x, 0) versus z line becomes 
universal, one is automatically finding the value of x be- 
yond which a saturation should be reached in the 8 ver- 
sus x plot. It is also a straightforward conclusion of the 
data in Fig. ^ that BSCO has a behavior very similar 
to YBCO although there is no direct evidence for in- 
commensurability or stripes in BSCO. It is interesting 
to notice that T c (x, z)/T c (x, 0) when plotted versus the 
in-plane 2D residual resistance, Pq D , exhibits universal 
behavior in the under-doped but not in the over-doped 
regime ^]. However, as a function of z, as we have shown 
here, the universal behavior is found in the optimally- 
and over-doped regimes. This fact is due to the linear 
dependence of Pq D with z || . 

We finally would like to point out that the predic- 
tion that z c tx T^{x,z — 0) is highly unusual. In 
BCS-like models the Abrikosov-Gor'kov |36) approach 
predicts that T c is suppressed when the scattering rate 
1 /t s oc z becomes of the order of the superconducting gap 
oc T c (x, 0) and therefore one would expect z c oc T c (x, z = 
0) in disagreement with the data. 

In conclusion, we presented a "geometrical" model for 
the suppression of T c (x,z) by planar impurities within 
the stripe model. It is based on the assumption that Zn 
impurities work as local pinning centers and suppress the 
transverse kinetic energy density of the stripes in their 
immediate vicinity. Our model leads to a linear T c sup- 
pression with doping, with a slope inversely proportional 
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to the critical impurity concentration z c , as given by Eq. 
(0). Moreover, we find that z c exhibits a non-trivial be- 
havior in under-doped compounds, z c (x) oc T£(x), and a 
universal one, z c =const. at optimally- and over-doped 
compositions. The obtained results agree very well with 
the experimental findings and allow to establish a con- 
nection between neutron scattering and magnetization 
measurements. 
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FIG. 4. NMR and 

(j,SH data for ytrium-and bismuth-based compounds. Red 
points correspond to YBa2Cu307-i @ co-doped with Li 
(squares) and Zn (circles), whereas open symbols correspond 
to Bi2Sr2Cai_ a Yj,(Cui_ z Co z )208+2, Q. (a) optimally- and 
over-doped compounds; (b) under-doped materials. 
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